Accepted by the Astro-physical Journal 



Four Highly Dispersed Millisecond Pulsars Discovered in the Arecibo PALFA 

Galactic Plane Survey 

F. Crawford 1 , K. Stovall 2 ' 3 , A. G. Lyne 4 , B. W. Stappers 4 , D. J. Nice 5 , I. H. Stairs 6 , P. 
Lazarus 7 ' 8 , J. W. T. Hessels 9 ' 10 , P. C. C. Freire 8 , B. Allen 11 ' 12 , N. D. R. Bhat 13 , S. Bogdanov 14 , 
A. Brazier 15 , F. Camilo 14 , D. J. Champion 8 , S. Chatterjee 15 , I. Cognard 16 , J. M. Cordes 15 , J. S. 

Deneva 17 , G. Desvignes 8,16 , F. A. Jenet 2 , V. M. Kaspi 7 , B. Knispel 11 , M. Kramer 8 , J. van 
Leeuwen 9 ' 10 , D. R. Lorimer 18 , R. Lynch 7 , M. A. McLaughlin 18 , S. M. Ransom 19 , P. Scholz 7 , X. 

Siemens 12 , A. Venkataraman 17 

ABSTRACT 

department of Physics and Astronomy, Franklin and Marshall College, P.O. Box 3003, Lancaster, PA 17604, 
USA; email: fcrawfor@fandm.edu 

2 Center for Gravitational Wave Astronomy, University of Texas at Brownsville, Brownsville, TX 78520, USA 

3 Department of Physics and Astronomy, University of Texas at San Antonio, San Antonio, TX 78249, USA 

4 Jodrell Bank Centre for Astrophysics, University of Manchester, Manchester M13 9PL, UK 

5 Department of Physics, Lafayette College, Easton, PA 18042, USA 

department of Physics and Astronomy, University of British Columbia, 6224 Agricultural Road, Vancouver, BC 
V6T 1Z1, Canada 

department of Physics, McGill University, 3600 University Street, Montreal, QC H3A 2T8, Canada 
8 Max-Planck-Institut fur Radioastronomie, auf dem Huegel 69, 53121 Bonn, Germany 

9 ASTRON, the Netherlands Institute for Radio Astronomy, Postbus 2, 7990 AA, Dwingeloo, The Netherlands 

10 Astronomical Institute "Anton Pannekoek," University of Amsterdam, Science Park 904, 1098 XH Amsterdam, 
The Netherlands 

11 Albert-Einstein-Institut, Max-Planck-Institut fur Gravitationsphysik, D-30167 Hannover, Germany, and Institut 
fur Gravitationsphysik, Leibniz Universitat Hannover, D-30167 Hannover, Germany 

12 Physics Department, University of Wisconsin at Milwaukee, Milwaukee, WI 53211, USA 

13 Swinburne University, Center for Astrophysics and Supercomputing, Hawthorn, Victoria 3122, Australia 

14 Columbia Astrophysics Laboratory, Columbia University, New York, NY 10027, USA 

15 Astronomy Department, Cornell University, Ithaca, NY 14853, USA 

1B Laboratoire de Physique et Chimie de l'Environnement et de l'Espace, LPC2E, CNRS et Universite d'Orleans, 
and Station de radioastronomie de Nangay, Observatoire de Paris, France 

17 Arecibo Observatory, HC3 Box 53995, Arecibo, PR 00612, USA 

^Department of Physics, West Virginia University, Morgantown, WV 26506, USA 

19 National Radio Astronomy Observatory, 520 Edgemont Rd., Charlottesville, VA 22903, USA 



- 2 - 



We present the discovery and phase-coherent timing of four highly dispersed millisec- 
ond pulsars (MSPs) from the Arecibo PALFA Galactic plane survey: PSRs J1844+0115, 
J1850+0124, J1900+0308, and J1944+2236. Three of the four pulsars are in binary sys- 
tems with low-mass companions, which are most likely white dwarfs, and which have 
orbital periods on the order of days. The fourth pulsar is isolated. All four pulsars 
have large dispersion measures (DM > 100 pc cm -3 ), are distant (> 3.4 kpc), faint 
at 1.4 GHz (< 0.2 mJy), and are fully recycled (with spin periods P between 3.5 and 
4.9 ms). The three binaries also have very small orbital eccentricities, as expected for 
tidally circularized, fully recycled systems with low-mass companions. These four pul- 
sars have DM/P ratios that are among the highest values for field MSPs in the Galaxy. 
These discoveries bring the total number of confirmed MSPs from the PALFA survey to 
fifteen. The discovery of these MSPs illustrates the power of PALFA for finding weak, 
distant MSPs at low-Galactic latitudes. This is important for accurate estimates of the 
Galactic MSP population and for the number of MSPs that the Square Kilometer Array 
can be expected to detect. 

Subject headings: pulsars: general — pulsars: individual (PSR J1844+0115, PSR 
J1850+0124, PSR J1900+0308, PSR J1944+2236) — surveys 



1. Introduction 



The PALFA survey is an ongoing, large-scale pulsar survey of the Galactic plane that uses the 
Arecibo 305-m radio telesc ope and the Arecibo L-Band Feed Array (ALFA) 7-beam multi-beam 
receiver (ICordes et al.ll2006l ). It is one of Arecibo's key science projects, and it will ultimately cover 
the entire Arecibo-visible sky within 5° of the Galactic plane (longitudes of 32° < I < 77° and 
!68° < I < 214°). PALFA observes at relatively high observing frequencies (1220 - 1520MHz) 
in order to mitigate the deleterious effects that interstellar dispersion and scattering have on the 
detection of distant pulsars at low Galactic latitudes. In this sense, the PALFA survey is similar 
to the highly successful Parkes Multibeam Pulsar Survey (PMPS; Manchester et al. 2001), but 
with increased tim e and frequency resolution, such as those of the current Parkes HTRU survey 
((Keith et alJhoid ). 



Compared to past Arecibo surveys, PALFA explores a far larger spatial volume due to its high 
time and frequency resolution, enabling discovery of faint, highly dispersed millisecond pulsars 
(MSPsQ in the Galactic plane, which tend to be at larger distances. Since the majority of MSPs 
have binary companions (~ 80%, Lorimer 2008), they are often interesting test cases for studies of 
exotic binary stellar evolution (e.g., Archibald et al. 2009; Freire et al. 2011). The unsurpassed 



20 An MSP typically refers to a (partially or fully) recycled pulsar having a small surface magnetic field strength 
(B < 10 10 G) and a large characteristic age (> 10 9 yr). 
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sensitivity of the Arecibo telescope is highly advantageous for detecting binary MSPs because the 
PALFA pointing dwell times are only a small fraction (less than 10%) of the orbital periods of 
all known binary radio pulsars. In this re gime, linear acceleration se arches are highly effective at 
recovering Doppler-smeared periodicities (j Johnston fc Kulkarnilll99ll ). 



The discovery of PSR J1903+0327 by IChampion et al.l (120081 ) is an excellent illustration of 
the PALFA survey's sensitivity to highly dispersed pulsars. PSR J1903+0327 has a spin period 
of only 2.15 ms and has the highest dispersion measure (DM) of all known completely recycled 
Galactic MSPs, 297 pc cm -3 . It occupies a region of DM-period phase-space that has previously 
been unexplored for pulsars in the Galactic field. For example, of the 90 recycled Galactic field 
radio pulsars currently listed in the ATNF catalog (apart from the PALFA discoveries jHl with spin 
periods P < 25 ms, only 9 have DM > 100 pc cm~ 3 , and only 14 have DM-derived distances listed 
in the catalog that are greater than 3 kpc. 

One of the main motivations for finding distant, highly dispersed MSPs is to determine a more 
complete census of the Galactic MSP population, which is currently biased by the large number of 
nearby (< 2 kpc) sources - especially given the recent discoveries of more than 40 generally nearby 
MSPs through targeted searches of Fermi gamma-ray sources (e.g., Ransom et al. 2011). MSPs are 
the longest-living active manifestations of neutron stars, with active lives hundreds to thousands of 
times longer than those of normal pulsars, magnetars or accreting neutron stars. Hence they give 
valuable insight into the Galactic neutron star population and binary stellar evolution in particular. 
The planned Square Kilometer Array (SKA; e.g., Carilli & Rawlings 2004) should be able to detect 
a large fraction of the MSPs in the Galaxy, and finding high-DM MSPs beforehand will help tell us 
how many to expect. Furthermore, MSPs are excellent probes of the interstellar medium, and the 
discovery of more distant MSPs whose signals likely pass through multiple scattering screens opens 
new possibilities (and challenges) in this area. For instance, scattering measurements of large-DM 
MSPs can be used to compare the observed effects of scattering on timing behavior with predictions 
of those effects. 

MSPs can also be used as high-precision astronomical clocks. Stable MSPs form the basis of 
the pulsar timing efforts of the International Pulsar Timing Array co nsortium to detect long-period 



gravi tational radiation from observations of pulsar timing residuals (IJenet et al 



2005; 



Hobbs et al 



201u ) . The re cent precision mass m easurements of PSR J1903+0327 (jFreire et al.ll201ll ) and PSR 
J1614— 2230 (IDemorest et al.l |2010| ) indicate that MSPs can have masses well above the Chan- 
drasekhar mass, and mass constraints from new MSPs will continue to map out the MSP mass 
distribution. The PA LFA pulsar PSR J1 949+3106, which may also have a higher mass than the 
Chandrasekhar mass (IDeneva et al.ll2012l ). is an example where future mass measurements may 
prove important. Although high-DM MSPs have timing precision problems associated with in- 
terstellar scattering in addition to their being generally radio- faint, future observations with next 



1 We also exclude the high-DM MSPs in globular clusters, because these are found by targeted searches conducted 
with very long integration times and at even higher observing frequencies (generally ~ 2 GHz). 
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generation instruments like the Square Kilometer Array (SKA) may be able to mitigate these fac- 
tors. Precision measurement of NS masses (e.g., Demorest et al. 2010) and the measurement of 
ultra-high spin rates (e.g., Lattimer & Prakash 2007; Hessels et al. 2006) can also rule out some 
high-density NS equations of state. 

In this paper we report the discovery and follow-up timing of four highly dispersed MSPs from 
the PALFA survey: PSRs J1844+0115, J1850+0124, J1900+0308, and J1944+2236. These MSPs 
all have DMs that are in the top 5% of non-PALFA field radio MSPs. Three of the four pulsars are in 
binary systems with low-mass companions. These di scoveries bring the nu mber of PALFA recycled 
pulsar discov eries to 15, including PSR J1903+0327 (jchampion et al-lhoO&h . PSRs J1949+3106 and 
J1955+2427 dpeneva et alJ boid). and the partially recycled PSRs J2007+2722 and J1952+2630 



(|Knispel et al.1 1201(1 



2011 



). The six additional recycled pulsars that have been confirmed in the 
survey need to be timed further to establish their rotational and orbital characteristics. These 
six pulsars will be published in forthcoming papers and are not considered further here. In §2 
we describe the discovery and follow-up observations of the four pulsars, and we present their 
phase-coherent timing solutions. In §3 we discuss our results, and in §4 we present our conclusions. 



2. Discovery, Timing, and Polarimetry Observations 



2.1. Discovery 

Until February 2009, survey data for the PALFA survey wer e recorded with the Wideband 



Arecibo Pulsar Processor (WAPP) auto-correlation spectrometers (jDowd et al.ll2000l ; ICordes et al 
2006 )1^1 The WAPP backends had 3-level sampling and provided 100 MHz of bandwidth for each 



of the 7 ALFA beams, centered on a sky frequency of 1.4 GHz and split into 256 lags (channels). 
These lags were sampled at 64 /us, and each survey pointing was observed for 268 s. The four 
pulsars presented here were discovered in data taken with this setup. 

PALFA survey data are archived at the Cornell University Center for Advanced Computing, 
where they are processed locally using custom search software and are also staged for transport to 
other PALFA consortium sites. The data are processed at these sites with dedicated computer clus 



ters t hat use a search pipeline based on the PRESTO software suite (IRansom 



2002 



23 



2001 



Ransom et al 



Another analysis pipeline is also used which employs the huge volunteer computing re- 
sources provided by the "Einstein@Home" project@ to search the survey data for binary systems 



22 Since 2009, the survey has been recording data using the Jeff Mock spectrometers, which are polyphase filterbanks 
that provide better radio frequency interference rejection as well as increased bandwidth compared to the earlier 
WAPP survey data. 

23 http: / /www.cv. nrao.edu/~sransom/presto 

24 r 



http://einstein.phys.uwm.edu 
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with orbital periods as short as 11 minutes ([Knispel et al.ll201ll ). Results from the PRESTO pro- 



cessing (including candidate plots and associated information) are uploaded to a collaborative web 
portal (www.cyberska.org), where interactive applications are used to visually inspect and rate 
selected candidates. 

The four MSPs presented here were all discovered with the PRESTO search pipeline in survey 
pointings taken at different epochs (see Table [1] for the discovery observation MJDs). Figure Q] 
shows the integrated pulse profiles for the four MSPs obtained by phase-aligning and adding the 
folded profiles from Arecibo 1.4 GHz observations. All of the integrated profiles in Fig. [1] have 
moderate widths (~ 15-30% of the pulse period). In two cases (PSRs J1844+0115 and J1850+0124), 
the profiles are clearly single-peaked. For PSR J1900+0308, there is some indication of a weaker 
leading component (and possibly also a trailing component). PSR J1944+2236 may have a hint of 
a secondary peak close to the center of the profile. None of the profiles shows any clear eviden ce of 



scattering, which is not surprising: according to the NE2001 model of ICordes &: Lazid (|2002l ). the 



pulse scattering t ime for these pulsars is expected to be small in all cases at 1.4 GHz (< 0.04 ms). 



Bhat et al.l ((2004) have noted that that in some cases pulse scattering times can be underestimated 
by the NE2001 model by up to an order of magnitude (see their Fig. 6). However, even if this were 
the case here, the scattering time would still be < 10% of the pulse period in all cases. At the large 
DMs of these pulsars, any flux variability would likely be from intrinsic mechanisms rather than 
scintillation. We see no significant variability or other intermittency (e.g., eclipsing effects) in the 
observations. 



2.2. Timing 

Soon after the discovery of these pulsars, follow-up timing observations began with the Arecibo 
telescope and the Lovell telescope at the Jodrell Bank observatory. The Arecibo observations used 
several systems: the ALFA receiver and single-pixel L-wide receivers, with data recorded with either 
the WAPPs or the Mock spectrometer. 

For the Jodrell Bank observations, dual-polarization cryogenic receivers covered 384 MHz of 
bandwidth centered at 1520 MHz. The bandwidth was split into 0.5 MHz channels in a digital 
filterbank, and each channel was sampled to provide 1024 samples per pulsar period. The effective 
sampling time for our four pulsars was therefore a few /is and depended on the pulsar period. Each 
pulsar was observed for between 2400 and 7200 s per timing session, depending on the pulsar. The 
sampled data in each channel were folded at the topocentric pulsar period in each case (this was 
predicted by an ephemeris), and the channels were subsequently dedispersed at the ephemeris DM 
and summed. 

For the Arecibo observations, the WAPP systems were used which covered two adjacent 50 
MHz bands. Each WAPP had 512 lags and a sampling time of 64 /xs. The pulsars were typically 
observed for 600 s in each observation. For PSR J 1900+0308, some observations used the Mock 
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spectrometers (see, e.g., Deneva et al. 2012). 

The folded profiles from the Arecibo observations were phase-aligned and summed to produce 
a high signal-to-noise pulse template for each pulsar (see Fig. [[]). Each observation was dedispersed 
and folded with the latest timing ephemeris, and folded profiles from each observation were cross- 
correlated in the Fourier domain with the profile template to obtain times-of-arrival (TOAs) with 
uncertainties. A separate pulse template was produced from the Jodrell Bank observations and 
this was used to similarly produce the Jodrell Bank TOAs. Table [1] shows the number of TOAs 
generated from Arecibo and Jodrell Bank with their typical rms values. 



We fit the resulting TOAs using the TEMPO software packager 5 ! and standard pulsar timing 
procedures (see, e.g., Freire et al. 2011 for details). Initial phase-connections for the pulsars were 
obtained from the Jodrell Bank observations alone, and these were subsequently supplemented 
with Arecibo TOAs at a variety of frequencies. The Arecibo TOAs allowed us to constrain the 
DM for each pulsar. All three binaries have small values for the product of orbital ecc entricity and 



projec ted semi-major axis (see Table [T]), so we used the ELL1 binary orbital model of lLange et al 



(2001) for the timing solutions. 

Figures [2] and [3] show the timing residuals for the four MSPs as a function of both date and 
orbital phase (except for PSR J1944+2236, which is isolated). As commonly happens in MSP 
timing, we found the formal TOA uncertainties in some data sets to be underestimated by up to 
a factor of 2.0. We scaled the TOA uncertainties by a common factor for each observational setup 
(telescope and backend combination) in order to make the reduced x 2 °f each TOA subset equal 
to unity, and we used the scaled uncertainties to weight the data in the timing fit. Increasing all 
uncertainties by a factor until the normalized x 2 is one is a conservative step, which produces less 
precise timing parameters, but ones that are justified by the actual rms of the TOAs available. 
The scaling factors are listed in Table [TJ There are no obvious systematic trends seen in the timing 
residuals, and the typical residual rms values from the timing solutions are less than 2% of the 
pulse period in each case (Table [T]) . These results indicate that the model is adequately describing 
the TOAs, with no significant unmodeled effects present. 

The full timing solutions for the four MSPs are presented in Table [IJ with the listed uncer- 
tainties representing twice the formal uncertainties produced by TEMPO. Apart from the directly 
measured astrometric, spin, and orbital parameters, derived physical parameters are also included. 
The physical parameters assume a pure magnetic dipole spin-down. Estimated distances and lu- 
minosi ties for the puls ars are also calculated and presented in the table. Note that the Shklovskii 



effect (I S hklovskiil 11970 ) and accelerations from the Galactic potential are not accounted for in the 
measured values of P. Until proper motions can be measured (see discussion below), we cannot 
correct for this effect. 



25 http://tempo.s ourceforRe.net. Note that we also obtained consistent results with the TEMP02 package 
( Hobbs et ai]|2006l ). 
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2.3. Polarimetry 

All four pulsars were also observed in a set of 1.4 GHz Arecibo polarimetry observations in 
order to measure polarization characteristics in the pulse profiles and to obtain rotation measures 
(RMs) and calibrated flux density estimates for the pulsars. Each pulsar was o bserved for betw een 



5 and 15 minutes at a center frequency of 1412 MHz using the ASP backend (IDemorestl 120071 ) . A 
total bandwidth of 24 MHz was split into 6 channels, each of 4 MHz, and the data were coherently 
dedisp ersed and fo l ded d uring the observation. The data were then processed using PSRCHIVE 



tools (jHotan et al.1 |2004j) . The resulting calibrated data files contained full Stokes parameters 
for each pulse profile bin. The Stokes parameters were converted to total intensity, linearly and 
circularly polarized intensity, and position angle (PA) values. 

We attempted to search for the RM for each pulsar in order to correct the data for Faraday 
rotation. We used trial RM values to see which trial produced the maximum linear polarization 
when the frequency channels were summed. This would indicate that the correct RM had been 
found. We were unable to find a reliable RM estimate in any of the four cases. Figure [4] shows 
the polarization profiles for the four pulsars without any Faraday rotation correction. None of 
the pulsars shows clearly significant polarization, with the possible exception of PSR J1900+0308, 
which may have an excess of right-handed circular polarization across the on-pulse bins. None of 
the pulsars has a linearly polarized intensity that is significant enough for reliable PA measurements 
across the on-pulse bins (in fact, only PSR J1844+0115 has any measurable PAs at all). Faraday 
smearing across the band is unlikely to be a significant factor in the small measured polarized 
signal. For RMs as high as a few hundred rad m -2 , the uncorrected Faraday rotation would reduce 
the linear polarization by only a few percent. 

The calibrated files were also used to determine 1.4 GHz flux densities for the pulsars (Table 
[1]). In all cases except PSR J1944+2236 we obtained flux density estimates, but only one digit of 
precision is quoted owing to the significant uncertainties in these estimates. For PSR J1944+2236, 
only an upper limit of 0.1 mJy is quoted for the flux density since a value was not measured. 
This was chosen as a reasonable upper limit since this is the value measured for the next weakest 
pulsar, and if PSR J1944+2236 had this flux density it would have likely been measurable. Longer 
observations in the future using the ASP or different observations using a system with a wider 
bandwidth may be useful in obtaining RM measurements and a clearer indication of the polarization 
characteristics of these MSPs. 



3. Discussion 

The inferred spin-down luminosities, surface magnetic field strengths, and characteristic ages 
of these four pulsars (Table [T]) are typical of the values of most fully recycled pulsars. Likewise, 
the orbital properties of the three binaries presented here are typical, and they have low-mass 
companions (M c ~ 0.2M Q ) and low eccentricities (e<3x 10 -4 in all cases), which suggests that 
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the systems are fully recycled with white dwarf companions. 

The measured eccentricities of the three binaries can be compared to the relat i onship be- 
tween the binary orb i tal pe riod, Pb, and eccentricity, e, that was outlined by iPhinnevI (|1992l ) and 
Phinney Kulkarnil (| 19941 ) for stable mass transfer from a Roche-lobe filling red giant (see also 
Camilo et al. 2001 and Lorimer 2008). PSRs J1844+0115 and J1850+0124 have eccentricities that 
easily fall within the predicted range of the model for their orbital periods. PSR J1900+0308 has 
an eccentricity that is small compared to the predicted range but is still consistent with the model. 

PSRs J1850+0124 and J1900+0308 may be useful f or tests of the Strong Equivalence Principle 
(SEP) and similar dev i ations from general r elativity (jDamour &; Schaeferl Il99ll ; iBell &: Damour 
1990 : Istairs et al.1 12003 : faonzalez etahlhoill ). The figure of merit for SEP tests is Pjf/e, which 
for PSR J1850+0124 is large, > 10 8 day 2 . T his places it in the sam e range as PSRs J1711-4322, 
J1933— 6211, and J1853+1303 (see Table 1 of lGonzalez et ajj (|201ll )). which were among the best 
pulsars used to constrain violation of the SEP a nd the strong-field version of the Parameterized Post- 
Newt onian (PPN) parameter 03 in that paper (|Will Nordtvedtll 1972l ; iDamour &: Esposito-Farese 
1992s! ). 

PSR J1900+0308 has a very large lower limit f or P^ 3 /e, making it potentially a key pulsar 



for co nstraining the (strong-field) PPN parameter d\ (IDamour &: Esposito-Faresdll992bl : IBell et al 



19961 ). Historically, the pulsars used for this test have had orbital periods of just a few days, since 
these tend also to have the lowest orbital eccentricities. Unfortunately, the short orbital periods 
of the relevant pulsars make it hard to define a population of systems thought to have followed 
the same general evolutionary path; in particular, pulsars with white-dwarf comp anions and orbital 



periods under about 4 days cannot safely be assumed to follow the Pb-m,2 relation (jRappaport et al 



1995k iTauris k. Savoniid Il999l ; iThorsett fc Chakrabartvl Il999i ) . Since a well-defined population for 
which one can make reasonable evolution-ba sed guesses at the pulsar and companion masses is 
desirable to mitigate against selection effects (jWexl l2000l ) . it has been difficult to define a suitable 
set of pulsars for use in the d?i test. The existence of PSR J1900+0308, assuming it does follow the 
Pb-mi relation, opens the possibility of using this class of pulsar for the d\ test. This will best be 
done once the proper motion and eccentricity are well-measured for this system. 

We are not currently able to measure proper motions for these pulsars. Making the assumption 
of a typical transverse speed of ~ 100 km s _1 and a representative distance of 5 kpc, the proper 
motion would be ~ 4.2 mas yr -1 . This would introduce a timing residual of ~ 5 fis after one year 
(jShklovskiil Il970l ) . and the timing spans for these pulsars are relatively short (~ 2 yr for three of 



the four pulsars; see Table H]). The relatively large rms values of the residuals (of order tens of fis) 
preclude these objects from being useful in pulsar timing arrays in the search for gravitational waves. 
Given their large distances, their proper motions are likely to be small and difficult to measure given 
their timing precision, making it hard to accurately determine the kinematic contributions to P 
and Pb. 



We performed a test with the timing residuals in order to assess the role of red noise (if 
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any) in the timing behavior of the four pulsars. The test we employed was the zero-crossing 
test, which measured the number of times the residuals change sign (cross zero) as a function of 
time (e.g., Deneva et al. 2012). The validity of this test is independent of the spacing of the 
individual timing points and corresponding residuals. The expected number of zero crossings of N 
residuals if only white noise were present is (Zyy) = (N — l)/2, with an uncertainty in the mean 
of (Jz w = \j (N — l)/2. The presence of timing noise (red noise) will produce fewer crossings than 
expected. In all four cases, the observed number of residual crossings fell within the expected range 
(Z\y) ± crz w , indicating that red noise is not significant in the timing behavior of these pulsars at 
the current (rather poor) level of precision. We have also char acterized the timing noise in these 4 



MSPs using a relative timing noise parameter, outlined by IShannon Corded (|2010l ). £ is the 
ratio of the observed timing residual value to a value expected from only timing noise. This latter 
parameter is determined using a scaling law (see their Eq. 7). The parameters used in the scaling 
were derived from measurements of canonical and millisecond pulsars (see Table 1 of Shannon & 
Cordes 2010). For all four MSPs, we found that this ratio was ( > 1. This indicates that white 
noise is dominating the residuals in all cases, and it is consistent with the results of the zero-crossing 
test above. 

The relatively large distances to these pulsars are restrictive in terms of their potential for 
multi-wavelength follow-up, and, unsurprisingly, the SIMBAD databas^fl indicates that there are 
no optical or IR counterparts listed at the positions of two of the three binary pulsars. There 
is an unidentified IR source, IRAS 18421+0112, located ~ 22" from PSR J1844+0115 which has 
a position uncertainty of 17" along the major axis of its position ellipse, but there is nothing 
listed in the 2MASS catalog at this location. It seems unlikely that this IRAS source is the 
binary counterpart of PSR J1844+0115, since the pulsar shows no evidence of eclipses or other 
variability /intermittency in the timing residuals that we might expect from a non-degenerate, 
extended companion (see Fig. [3]), and we see none of the timing jitter that is clearly present in the 
two known cases of these kinds of Galactic MSP systems that have been published to date (PSR 
J1023+0038, Archibald et al. 2009; and PSR J1723-2837, Crawford et al. 2010). 

There are also no known X-ray or 7-ray counterparts in the HEASARC catalog at the four 
pulsar positions. This is not surprising since E < 10 34 erg s _1 and for this spin-down luminosity 
range, the X-ray emission would be quite faint. To search for 7-ray pulsations from the four MSPs, 
we retrieved Fermi LAT data from the start of the mission up to 2012 February 6 using 1° extraction 
radii centered on the pulsar positions. The event lists were filtered using the recommended cuts in 
maximum zenith angle (100°), event class (2), and photon energy (> 100 MeV). The 7-ray events 
were folded using the Fermi plug-in for TEMPO_E_] and the radio ephemerides from Table 1. In all 
instances, no statistically significant pulsations were detected. Repeating the analysis with a >300 
MeV energy cut and extraction radii in the range 0.5° — 1° also yielded no detections. 



'http://simbad.u-strasbg.fr/simbad/ 



See |http7/7 fermi.gsfc.nasa.gov/ssc/data/analysis/user/Fermi_plug_doc.pdf 
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All four of these MSPs share the quality that they have very large DMs and DM-inferred 
distances. This is also a feature of the other MSPs that have been discovered by the PALFA 
survey. Previous large-scale surveys had poorer sensitivity to these kinds of MSPs owing to their 
lower observing frequencies and inadequate observing instrumentation (wider frequency channels 
and insufficient sampling rates). Figure [5] illustrates this with a plot of DM vs. spin period for 
all 90 non-P ALFA Galactic held rad io pulsars in the ATNF pulsar catalog having P < 25 ms and 
P < 10~ 17 f Manchester et alJ 12005 1 28 L This group of recycled pulsars does not include globular 



cluster (GC), radio quiet, or young pulsars. Also shown ar e the four MSPs described here and 



the five other PALFA MSPs th at have been reported to date ([Champion et al.ll2008l ; iKnispel et al 
2010. 



2011 



Deneva et al.ll2012l ). It is clear that the PALFA MSPs occupy an area of the DM-period 



phase-space where very few MSPs are currently known. 

Figures [6] and [7J also illustrate the power of the PALFA survey for finding fast, distant MSPs. 
Both figures include the same set of data as shown in Fig. [5j Fig. [6] shows a histogram of the 
ratio DM/P (excluding PSR J1903+0327, since it falls well beyond the plot limits). The subset 
of PALFA MSPs are indicated with shaded parts of the histogram. The positions of the 4 MSPs 
presented here are indicated with arrows, all of which fall near the edge of the histogram. This 
plot demonstrates how we are exploring a larger DM/P parameter space with this survey. Fig. [7J 
shows a Galactic projection plot of the same set of pulsars. The nominal locations of the PALFA 
MSPs in the Galactic plane are indicated by stars. Exploring the high DM/P parameter space is 
necessary for a full Galactic census of MSPs. 

To demonstrate the potential of the PALFA survey in such a population analysis, we have 
carried out a preliminary investigation to model DM/P for Galactic MSPs. We used the freely 
available PSRPOP software packagj^l to carry out this work, and we generated simple "snapshot" 
models of the population which are normalized to reproduce the yield of 20 MSPs detected in the 
PM PS (e.g. L o rimer et al. 2006). The simulation procedure is described in that paper and also 



m 



Smits et al.1 (|2009i ). In brief, we generated four different models (labeled A through D) of the 



present day MSP population which matched the number of MSPs detected by the PMPS, but 
produced a variety of different DM/P values and predicted yields for the PALFA survey. For the 
purposes of this work, we focused on changing assumptions about the spatial distribution of the 
underlying MSP population. 

Model A assumes a Gaussian radial density profile with a standard deviation &r = 6.5 kpc, an 
exponential scale height with a mean of h z = 500 pc, the log-normal 1400-MHz pulsar luminosity 
function foun d by iFaucher-Giguere fc Kaspil (J2006J), and the MSP period distribution and beaming 
model used in lSmits et al. ( 2009 ). Model B is similar to Model A with the exception that the radial 
density is uniform (i.e., a constant surface density) throughout a disk of radius 25 kpc. Models C 



28 http://www.atnf. csiro.au/research/pulsar/psrcat/ 
29 http://psrpop. sourceforge.net 
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and D both follow Model A, but with h z set to 250 and 750 pc, respectively. Further models should 
be explored in a subsequent analysis, but for the purposes of the present work, the above choices 
demonstrate the importance of DM/F as a diagnostic for population studies. 

The results of our simulations are shown in Table [2] where we tabulate, along with the observed 
sample, the numbers of detectable MSPs along with the median DM/P found in each simulation. By 
definition, due to the way we normalize the simulations, the number of detected MSPs in the PMPS, 
-ApmpSj is 20. For models A through D, we tabulated the predicted number of MSPs detectable in 
the PALFA survey, Apalfa- The median DM/P values for each simulation can be compared with 
those currently observed in the survey samples. As can be seen, there is a tremendous variation 
in the predicted values, and the median DM/P depends strongly on the choice of scale height and 
radial distribution. Model B appears to underestimate both the observed number of PALFA MSPs 
(assuming we find more as the survey proceeds) as well as the median DM/P value. Model C 
seems to overestimate the median DM/P value. Models A and D are most consistent with our 
observations to date, suggesting that a large scale height h z and a Gaussian radial profile density 
are favored. 

4. Conclusions 

With its high time and frequency resolution and relatively high observing frequency, the PALFA 
survey is sensitive to Galactic field MSPs at large distances and DMs. The survey has so far 
discovered a total of 15 MSPs with P < 25 ms, 14 of which have DM > 100 pc cm -3 . Four of 
these PALFA MSPs (PSRs J1844+0115, J1850+0124, J1900+0308, and J1944+2236) are presented 
here with phase-coherent timing solutions. As an ensemble, the PALFA MSP discoveries show the 
ability of PALFA to extend the volume of MSP discovery space to relatively high DMs. This is the 
first step toward a nearly complete census of the Galactic MSP population, which will be possible 
with the SKA. 
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dation for Innovation, and a Killam Research Fellowship. VMK holds the Lome Trottier Chair 
in Astrophysics and Cosmology and a Canadian Research Chair in Observational Astrophysics. 
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support of the Max Planck Society. PL was partly funded by a NSERC PGS scholarship and an 
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Fig. 1. — Integrated 1.4 GHz pulse profiles for four MSPs. The pulse profiles were constructed 
by phase aligning and adding between 30 and 50 minutes of Arecibo observations for each pulsar. 
A combination of WAPP and Mock data were used for the profiles. In all cases, 64 profile bins 
were used, except for PSR J1944+2236 where 128 bins were used. One full period is shown in each 
case. The vertical axis is arbitrary and has been scaled so that all profiles have the same maximum 
value. The horizontal bars indicate the amount of dispersion smearing within channels for Arecibo 
observations taken with the WAPPs. 
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Fig. 2.- 



Timing residuals for four MSPs as a function of date. Black points indicate Jodrell 
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Fig. 3. — Timing residuals for three binary MSPs as a function of orbital phase. PSR J1944+2236 
is isolated and is not shown. There are no systematic trends in the residuals. None of the three 
binaries show evidence of eclipsing effects which would be expected if the companions were extended, 
non-degenerate stars with large orbital inclination angles. This suggests that the companions 
are probably low-mass white dwarfs, though more massive companions with very small orbital 
inclinations cannot be ruled out. 
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Fig. 4. — Polarization profiles fo r the four MSPs taken at a center frequency of 1412 MHz with 
Arecibo using the ASP backend (IDemorestl 120071 ) . RMs for the pulsars could not be determined 
from these observations, so none of the profiles have been corrected for Faraday rotation. Each 
plot has 128 phase bins, except PSR J1944+2236 which has 64 bins. In the bottom part of each 
plot, the black, red, and blue lines correspond to total, linearly polarized, and circularly polarized 
intensity, respectively. The top part of each plot shows the linear polarization PA for bins with 
S/N > 3. The low signal strength in the observations is likely to be responsible for the lack of any 
significant linear polarization observed and any measurable PAs across the profiles. 
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Fig. 5. — Dispersion measure vs. spin period for 90 Galactic field radio MSPs from the ATNF 
pulsar catalog (only pulsars with periods less than 25 ms and P < 10 -17 are plotted, and GC and 
radio-quiet pulsars are not included). These pulsars are plotted as dots. Also plotted are the four 
PALFA MSPs reported in this paper (crosses) and five o t her MSPs discovere d by PALFA (filled 
circles) ([Champion et al.1 bood : knispel et aD bold , boill : beneva et all boid l. All nine PALFA 



pulsars are labeled. It is clear that the PALFA survey is exploring new parameter space with 
the discovery of these pulsars. The DMs of these pulsars are among the highest that have been 
previously discovered for MSPs in the Galaxy. 
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Fig. 6. — Histogram of DM/P for 90 known field radio pulsars from the ATNF catalog and 8 PALFA 
MSPs (not including PSR J1903+0327) having P < 25 ms, P < 10~ 17 . GC and radio-quiet pulsars 
are not included here. PSR Jf903+0327 is not shown owing to its very large DM/P value of 138.4 
pc cm~ 3 ms -1 , which is far off the scale. The unshaded histogram shows all 98 pulsars, including 
the 8 PALFA pulsars. The shaded histogram shows the subset of 8 PALFA MSPs. The four MSPs 
presented in this paper are indicated by arrows and are labeled. It is clear that the PALFA survey 
is finding fast, distant pulsars that were undetectable in previous large-scale surveys. 
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Fig. 7. — Overhead Galactic projection plot showing the nominal locations of 90 known field radio 
MSPs in the ATNF catalog with P < 25 ms and P < 1 ~ 17 (blue dots) and 9 kno wn PALFA MSPs 



(red stars) (jChampion et al.1 120081 ; iKnispel et al 



2010 



2011 



Deneva et 12012 ) . including the 4 



described in this paper. The nominal distances are inferred from the DMs and are uncertain by 
~ 30%. The Sun is located at the center of the cluster of dots, and the Galactic center is at the 
origin. The black curves indicate locations of the Galactic spiral arms according to the NE2001 
model of Cordes & Lazio (2002). The cones indicated by the two pairs of straight lines indicate 

thp rpffinns nf thp dalnptip nlanp that arp visihlp frnm Arppihn Thp PAT.FA MSPs arp amnmr thp 
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Table 1. Timing Parameters for Four MSPs 



PSR 


J1844+0115 


J1850+0124 


J1900+0308 


J1944+2236 


Right ascension (J2000) 


18:44:40.5474(3) 


18:50:01.0139(8) 


19:00:50.5548(3) 


19:44:01.0707(11) 


Declination (J2000) 


+01:15:34.974(11) 


+01:24:34.61(3) 


+03:08:24.079(13) 


+22:36:22.62(2) 


Spin frequency, / (s — **) 


238.91757323111(15) 


280.9175173300(4) 


203.6975581461(2) 


276.3963275801(5) 


Frequency derivative, / (s 2 ) 


—6.12(12) x 10~ 16 


—8.6(4) x 10~ 16 


— 2.45(9) X 10~ 16 


— 5.7(4) x 10~ 16 


Dispersion measure, DM (pc cm~^) 


148.22(2) 


118.89(5) 


249.898(11) 


185.45(12) 


Orbital period, P(, (d) 


50.6458881(11) 


84.949858(4) 


12.47602144(10) 




Time of pcriastron passage, Tq (MJD) a 


55409.21(4) 


55241.3(3) 


_ 


_ 


Projected semi-major axis, x (s)^ 


14.173495(9) 


34.00102(2) 


6.716377(4) 




Longitude of periastron, uj (dcg) a 


189.6(3) 


98.8(11) 






Eccentricity, e 


2.578(11) X10~ 4 


6.90(12) X10~ 5 


< 3.3 X 10~ 6 




Time of ascending node, T asc (MJD) C 


55382.532417(5) 


55217.963824(11) 


55304.460562(3) 




ei — e sin uj c 


-4.32(13) X 10~ 5 


6.82(12) x 10~ 5 


< 3.0 X 10~ 6 


- 


62 — e cos uj c 


-2.542(11) x 10~ 4 


-1.05(13) X 10~ 5 


< 3.1 X 10~ 6 




Period, P (ms) 


4.185543936664(3) 


3.559763768043(5) 


4.909239016417(4) 


3.617993078111(6) 


Period derivative, P 


1.07(2) x 10~ 20 


1.09(5) X 10~ 20 


5.9(2) x 10 -21 


7.5(5) X 10 -21 


Mass function (Mq) 


0.001191860(2) 


0.005848349(13) 


0.002089949(4) 




Companion mass (Mq) 


> 0.14 


> 0.25 


> 0.17 




Galactic longitude, I (deg) 


33.28 


34.02 


36.79 


58.90 


Galactic latitude, b (deg) 


+2.08 


+0.96 


-0.66 


-0.66 


Surface magnetic field, B (G) e 


2.1 x 10 s 


2.0 x 10 s 


1.7 x 10 s 


1.7 x 10 s 


Spin-down luminosity, E (erg s 1 ) c 


5.8 X 10 33 


9.6 X 10 33 


2.0 X 10 33 


6.4 x 10 33 


Characteristic age, T c (Gyr) c 


6.2 


5.2 


13.0 


7.5 


Distance, d (kpc) 


3.9 


3.4 


5.8 


6.5 


Distance from Galactic plane, \z\, (kpc) g 


0.14 


0.06 


0.07 


0.08 


1400 MHz flux density, S1400 ( m Jy) 


~ 0.1 


~ 0.2 


~ 0.2 


5 0.1 


1400 MHz radio luminosity, -L1400 ( m Jy kpc 2 ) h 


~ 1.5 


~ 2.3 


~ 6.7 


< 4.2 


Discovery observation MJD 


53493 


54555 


53656 


53647 


Timing epoch (MJD) 


55383 


55383 


55255 


55415 


TOA range (MJD) 


55152-55914 


55153-55915 


54907-55927 


55222-55926 


Timing span (d) 


762 


762 


1020 


704 


Number of points in timing fit ( Jodrcll/Arccibo) 


116/7 


110/7 


94/136 


77/8 


Characteristic TOA residual (^is) (Jodrcll/Arccibo) 


40/10 


100/10 


80/10 


80/10 


TOA error scale factor (Jodrell/WAPP/Mock)' 


1.20/1.00/- 


2.03/1.00/- 


1.09/1.25/1.45 


1.87/-/1.85 


Weighted rms post-fit residual {p.s) 


25.4 


51.6 


15.8 


66.1 



Note. — Figures in parentheses arc uncertainties in the last digit quoted and arc twice the formal errors from the TEMPO timing solution. 

a The parameters Tq and u arc highly covariant in the timing solutions. Observers should use the following values of Tq and uj, respectively: 
55409.21301273 and 189.6504286222 for PSR J1844+0115; 55241.26871848 and 98.7613423916 for PSR J1850+0124. 

b a; — asini/c where a is the semi-major axis and i is the orbital inclination angle. 

c For all three binaries, the extremely small eccentricity (e <g; 1) introduces a covariancc term between the time of pcriastron passage Tq and the 
longitude of pcriastron u. In these cases the ELL1 binary model was used jLange et al.1 2001). where the time of ascending node T asc , defined as 
when w — 0, and — e sin uj and 62 — ecosu are fit instead. Tq, e, and uj are derived from this. 

^Assumes an inclination angle i — 90° and a pulsar mass of 1.35 solar masses. 

C B = 3.2 X 10 19 (PP) 1//2 ; E = 47T 2 IP/P 3 , with an assumed moment of inertia / — 10 45 g cm 2 ; r c = P/2P. Note that these parameters depend 
on P which may be affected by the Shklovskii effect. 

f From the NE2001 DM-distance model of lCordes & Lazid j2002D . 

s \z\ — d sin 1 61 
h -Ll400 — Si40Qd 2 

"TOA uncertainties from each instrumental setup were multiplied by this factor (EFAC in TEMPO) to correct for the generally underestimated 
uncertainties produced by TEMPO (sec the text for justification). 
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Table 2. MSP Population Simulation Results for Four Models 



Sample 


-/VpMPS a 


ApALFA a 


DM/P PM Ps b 


DM/PpALFA b 


Observed 


20 


> 15 


8 


33 


Model A 


20 


40 


14 


56 


Model B 


20 


15 


9 


7 


Model C 


20 


37 


11 


113 


Model D 


20 


42 


16 


27 



Note. — The four models and their parameters are described 
in the text in §3. 

a Number of detectable MSPs. 

b Median DM/P value of detectable MSPs. 



